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Abstract 

Prompted by an example arising in critical percolation, we study 
some reflected Brownian motions in symmetric planar domains and show 
that they are intertwined with one-dimensional diffusions. In the case 
of a wedge, the reflected Brownian motion is intertwined with the 3- 
dimensional Bessel process. This implies some simple hitting distributions 
and sheds some light on the formula proposed by Watts for double-crossing 
probabilities in critical percolation. 



1 Introduction and notations 

It has been pointed out in [WerOll ILSW02) that the hitting distribution of a 
certain reflected Brownian motion (RBM) in an equilateral triangle was uniform. 
More precisely, if the Brownian motion is started from one corner and reflected 
on the two adjacent sides with the oblique reflection angle 7r/6 away from the 
normal direction toward the opposite side, then the Brownian motion will hit 
this opposite side with uniform distribution. Their proof uses a discrete version 
of this result (reflected simple random walk) on a well-chosen triangular lattice, 
and an invariance principle (i.e. the reflected random walk converges to reflected 
Brownian motion). This result is then combined with a locality property in 
order to identify the law of the whole "hull" of this stopped reflected Brownian 
motion with that of a chordal SLE 6 process (that has the same uniform hitting 
distribution and has also a locality property). 

In |Smi01j . Smirnov proved that SLE 6 is the scaling limit of critical perco- 
lation cluster interfaces on the triangular lattice. The main step in this proof is 
the derivation of the fact that in the scaling limit, the above-mentioned hitting 
distribution (if one replaces the RBM by the percolation interface) becomes 
uniform. It makes an important use of some specific features of the triangular 
lattice. Note that in the case of critical percolation, the discrete hitting prob- 
abilities are not uniform, as opposed to those for the reflected simple random 
walk. 
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One of the motivations of the present paper is to see whether one can gen- 
eralize this hitting probability property of reflected Brownian motion to other 
symmetric shapes (for instance other wedges) and other reflection angles. It 
turns out that one can choose the reflection angles in such a way that hit- 
ting probabilities of segments orthogonal to the symmetry axis of the shape 
remain uniform. Furthermore, this fact is closely related to a relation between 
these reflected Brownian motions and one-dimensional diffusions. More pre- 
cisely, the projection of the RBM on the symmetry axis is a diffusion in its 
own filtration; for wedges, one recovers the 3-dimensional Bessel process. This 
provides a new example of non-trivial intertwining relations (see for instance 
|RogPit8T| ICPY9 8 ) . Our proofs rely also on reflected simple random walks 
and an invariance principle. 

We shall see that the formula proposed by Watts |Wa96| for double-crossing 
probabilities (i.e. simultaneous top-to-bottom and left-to-right crossing of a 
quadrilateral) for critical percolation and for which it is not clear that the SLE6 
approach will confirm it, is in fact satisfied (in some appropriate sense) by 
the reflected Brownian motion; this will follow from the study of time-reversed 
reflected Brownian motions. 

This paper is organized as follows. First, we define some reflected random 
walks in wedges and study their scaling limit. Then we prove that the limiting 
reflected Brownian motions are intertwined with the 3-dimensional Bessel pro- 
cess. Making use of this fact, we determine the time reversals of these RBMs, 
and discuss an analogue of Watts' formula. Finally, we generalize some of the re- 
sults to symmetric shapes ("vases"), essential replacing the 3-dimensional Bessel 
process by a general one-dimensional diffusion. 



2 Invariance principle 

Let a £ (0, ^). We will consider the wedge 

C a — {z G C : —a < arg z < a} C C 
and the rectangular lattice 

L a = cos al + i sin olL. 

We will study random walks on the graph T a = L a n C a . 

In this section, the angle a is fixed, so we will drop the subscript a. Consider 
the following random walk on the graph T = L n C: for an inner point x, the 
transition probability is: 

p = p(x, x + cos a) = p(x, x — cos a) = sm 2 a 
q = p(x, x + i sin a) = p{x, x — i sin a) = cos 2 a 

so that such a step has a zero mean, and its covariance matrix is a multiple 
of the identity. If x ^ is a boundary point with arg(i) = a, the transition 
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Figure 1: A wedge with its associated rectangular lattice 



probability is: 

p(x, x + cos a) = p(x, x + cos a + i sin a) 
p(x, x — i sin a) = p(x, x) = cos 2 - 

Notice that in this case, 

BXg{W{X 1 -x))=2a-^. 

On the other boundary line, the transition probabilities are defined symmetri- 
cally. 

Finally for the apex, set 

r = p(0, cos a) = p(0, cos a ± i sin a) = - — ^ > 0. 

o 

The exact value of this positive probability will not matter in the continuous 
limit. This random walk starting from x £ T will be denoted by (X%) n >o. We 
will also call (Y t x )t>o the Brownian motion in C starting from x and reflected on 
the boundary with angle —a (see |VarWil85j for a definition of this process; we 
use the same conventions for reflection angles, i.e. positive angles point towards 
the apex). 

The following result is classical. We include a proof for the sake of complete- 
ness (and in order to point out later how it can be generalized). 

Lemma 1. Let x be a point in C , and (x n ) a sequence in T such that \ < 
1/2 for all n £N. Then the following weak convergence of processes holds as n 
goes to infinity: 

* (Xt )t>0 

t>0 



sin 2 a 

~ 2 



[n^t/ sin 2 a cos 2 a] 
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Figure 2: Transition probabilities for an inner point and a boundary point 

Proof. Let (Y^) be a continuous C- valued process that interpolates linearly 
the process: 

-X?" 
n 

Choose a Holder exponent j3, < /? < |. Then it is standard to check that for 
any fixed T > 0, the following estimate on /3-H61dcr norms holds : 

lim hmsupP(||rW|| [0i T],/3 > K) = 

One proceeds exactly as in the proof of the a.s. /3-H61der continuity for Brownian 
paths (there is no particular problem due to the boundary here) . Hence we get a 
tight sequence in the Polish space of /3-H61der continuous mappings from [0, T] 
to C. Using Prohorov's theorem, one gets the existence of subsequential weak 
limits. Then a standard diagonal argument yields subsequential weak limits for 
the whole process in the Wiener space of C-valued processes. 

So we have to check that there is only one possible weak limit, namely the 
reflected Brownian motion (Y t ). To prove this, one can use the submartingale 
problem characterization of reflected BM ( | Var Wil85] . Theorem 2.1). More pre- 
cisely, the RBM is the only C-valued process starting from x such that for any 
/ € C"l (C) (twice differentiable with bounded derivatives) with positive deriva- 
tives on the boundary along the reflection direction, the real valued process: 

f{Y t )- l -j\f{Y s )ds 

is a submartingale. So let / be such a function; it is sufficient to check that for 
any < s < t: 

liminf E (B) (f(Y t ) - f(Y s ) - I j Af(Y u )du^j > 

where E(„) designates the expectation operator for the n-th approximate pro- 
cess. Consider the stopping time r = inft^odY^I > M) for some large number 
M. It is enough to check the submartingale inequalities up to time r. Notice 
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that the following Taylor expansion holds: 
(pf(x + cosa/n) + pf(x — cosa/ro) + qf(x + i sina/n) + qf(x — i sma/n)) 

= —7 sin 2 a cos 2 aAf(x) + o(^r) 

and that the error term is uniform in x e C U D(0, M) (using for instance 
Hadamard's lemma). Let Uk = k sin 2 a cos 2 a/n 2 , k € N, correspond to the 
discrete jump times, and let B be the boundary of C . Then: 

~ f(Y sAT )) = E^ n) I E /( r ^ +1 ) - 

\s/\r<u k <tAr 

%)( E ^cos 2 asin 2 aA/(y; 

\sAr<u k <tAr 



+ e(Y Uk )U Yuk eBj +o(l) 
= Ef n) (£ AT l -Af{Y u )du^j + o(l) + E ^U(LyMy) 



y&B 

The error term o(l) is simply a Ricmann sum error. Then there is a "boundary 
error term" . In the formula, L y designates the local time at y (number of 
discrete jumps to y) between times s A r and t A r, and e(y) is defined as (if 
arg(y) = a for instance): 

e(y) = - f(y + isma/n))+p(f(y + e la /n) - f(y - cos 07^)) 

The first order term is proportional to the derivative of / at y along the reflection 
direction, hence is positive by hypothesis. One may write: e(y) + 0(l/n 2 ) > 
for y G B, the error term 0(l/n 2 ) being uniform in y G C U D(0, M). So the 
only thing to check now is that the time spent on the boundary is negligible, 
more precisely: 

EEf n) (i,) = o(n 2 ) 

y€B 

Recall that we consider the local time between times s and t, and before exit 
time r; one may drop this last condition for the sake of simplicity. Justifying 
the very intuitive fact that the walk spends a negligible time on a negligible part 
of the state space is rather tedious. The discrete intertwining relations allow 
explicit computations that yield the result when the starting point is the apex 
(see below) . Then, using the Markov property, one sees it is enough to prove 
this estimate for zero. This concludes the proof. □ 

One may carry out the same proof with different transition probabilities on 
the boundary, corresponding to different reflection angles; to derive the final 
local time estimate in the general case, one may use the local Central Limit 
Theorem. 
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Figure 3: "stretched" triangular lattice 



Remark 1. Similar results hold for different lattices, in particular "stretched 
triangular lattices", i.e. images of the standard triangular lattice under an or- 
thogonal affinity. In the special case a = n/6, one may use the standard tri- 
angular lattice, with the usual transition probability for inner point (i.e., 1/6 
for each neighbour). In the case a = tt/4, these two approaches give the same 
lattice, namely the square lattice. 

Note that all these reflected Brownian motions have the same parameter 
a = — 1 in Varadhan- Williams conventions ( VarWil85 ) (what we call a would 
be denoted £/2 in IVarWilSflQ . 

3 Intertwining Relations 

First, let us recall the notion of intertwined Markovian semigroups (for some 
background on the subject, see |CPY98| ). Let (Pt,t > 0) and (Q t , t > 0) be two 
Markovian semigroups with respective state spaces (Sp,Sp) and (Sq,Sq), and 
A a Markov transition kernel from (Sq,Sq) to (Sp,Sp). The two semigroups 
are said to be intertwined by A if the following identity of Markov transition 
kernels from Sq to Sp holds : 

APt = QtA 

for all i > (t is either a discrete or a continuous time parameter). 

Suppose that there exists a measurable function <f> : Sp —> Sq such that 
A(y, </> -1 (2/)) = 1 for all y e Sq. Then, as shown in |RogPit81| , if (Z t ) is a 
Markov process with semigroup (Pt) and initial law A(yg, .) for some j/o £ Sq, 
and Y t = 4>(Z t ), then (Y t ) is a Markov process with semigroup (Q t ) and starting 
state yo- Moreover, for any fixed time T, the following filtering formula holds: 

¥(Z T e.\Y t ,t>0)=A(Y T ,.). 
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Under appropriate regularity conditions (which will be satisfied in all exam- 
ples that we shall consider here), this formula also holds for almost sure finite 
stopping times T (w.r. to the filtration of Y). In particular, for T = T y , the 
first hitting time of y by Y , if T y < 00 a.s., then the random variable Zt v has 
distribution A(y, .) and is independent of (Y t ). 

Consider as before the random walk (X n ) on the graph L. Let A be the 
following Markov transition kernel from N to A (we omit to mention that these 
discrete spaces are equipped with discrete a- algebras): 



1 x 

A(x, .) — — j - ^ ^ $x cos a.-\-ki sin a 



k— — x 



Let (P n ) be the semigroup of (X n ), and let (Q n ) be the semigroup of the random 
walk on N with transition probability: 



sin a 2i — l 



q(i,i-l)= 2 2i+1 
q(i, i) = cos 2 a 

g ( i|i + l) = fi !^a«±3 

for i > and q(0, 0) = 1 — 3r, q(0, 1) = 3r. Note that i 1— > ^qrj- is harmonic for 
this transition kernel, except at 0. Then the intertwining relation holds: 

AP„ = Q n A 

for all n > 0. 

Let (Y t ) be the reflected Brownian motion in the wedge C with reflection 
angle —a. Let (P*)t>o be its semigroup. Let (Qt)t>o be the generator of the 3- 
dimensional Bessel process (with values in K+), and A be the Markov transition 
kernel from R + to C such that A(x, .) is the uniform distribution on the segment 
[x — i tan a, x + i tan a] . 

In the scaling limit, the previous intertwining leads to: 

Proposition 1. Then P and Q are intertwined by A. In particular, for x > 0, 
Z t = VtY t A{x '- ] is a 3- dimensional Bessel process starting from x. 

Proof. It is obvious that the intertwining relation holds in the limit. The only 
thing to check is that the intertwined random walk on N converges to a 3- 
dimensional Bessel process. One may proceed as in the previous proof, after 
computing: 



sin 2 a { 2nx/ cos a + 3 
2 y 2nx/ cos a + 1 
2nx/ cos a — 1 



cos af{x)-\ — — -f(x + cosa/n) 

\2nx/ c' 



2nx/ cos a + 1 



f(x -~ cosa/n) I I — f(x) 



sin 2 a cos 2 a / 1 . . 1 



x f'(x) + ^f"(x))+o(l/n 2 ) 
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Here one recognizes the generator of the 3-dimensional Bessel process (the factor 

•22 

Sln a n ? s " corresponds to the time scaling). Such discrete approximations of 
the 3-dimensional Bessel process are classical. □ 

The Markov transition kernel A acts on bounded Borel functions in the 
following way: if / is a bounded Borel function on C, then A/ is a bounded 
Borel function on R+ such that: 

(Af)(x) = E(/(x + iZ7tana)) 

where U is a uniform random variable on [—1,1]. One may remark that these 
intertwining relations fit in the "filtering type framework" described in |CPY98j . 

As Jim Pitman pointed out to us, there is an analogy with the situation in the 
(2M — X) theorem: let B be a standard (real) Brownian motion starting from 0, 
and let M t be its supremum up to time t. Then the process (2M t — X t ) t > is a 
3-dimensional Bessel process (in its own filtration) (see |RogPit81|lRevYor91p . 
A proof of this fact involves an intertwining relation. More precisely, note 
X t = Mt — B t , Yt = 2M t — B t . Then, if (Gt)t>o designates the natural filtration 
of (Y t )t>o, the following relation holds for every Borel function / : K + — > M + : 

E(f(X t )\g t )= [ f(xY t )dx 
Jo 

The Markov kernel in this last situation: Af(y) = f(xy)dy is very similar to 
the one we described earlier. 

The Brownian motion in a plane strip {z : < 3z < 1} with normal 
reflection on the boundary may be seen as a degenerate limiting case when a \ 
0. Indeed, the reflection direction makes an angle it/2 — a with the boundary, 
and the two boundary half-lines make an angle 2a. This particular RBM may 
be represented as Z = X + if(Y), where X and Y are independent standard 
(one-dimensional) BMs, and / is the "seesaw" function: for all k E Z, 

f(x) =x-2k if 2k < x < 2k + 1 

f(x) =2k + 2-x if 2k + 1 < x < 2k + 2 

In this case, the intertwining relation is easily proved. Indeed, let <p be any 
bounded Borel function on the strip and t be a fixed positive time. We have to 
check that: 

E x (cb(X t +iU))= E(F. x+iU (<t>(Z t ))) 

where U is an independent uniform random variable. Since the real and the 
imaginary part of Z are independent, we have only to check that, for any 
bounded Borel function tp on [0, 1]: 

E( 9 (U)) = E(Eu(<p{f(y t )))) = E (v(f(U + Y t ))) 

which is readily seen (if V is any symmetric r.v. independent from U, f(U + V) 
is a uniform r.v.). 
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Remark 2. In the general case, the intertwining relation does not seem obvi- 
ous directly in the continuous setting, hence the use of discrete approximations. 
Also, the real part of the RBM started from a fixed, inner point (not the apex) 
does not appear to have similar properties; in this case, the relation between the 
filtration of the RBM and the filtration of its real part seems very intricate, and 
the real part is generically not Markovian in its own filtration. 



4 Time reversal 

In this section, we consider the time reversals of the reflected Brownian motions 
studied in the previous sections. The intertwining relation enables an explicit 
computation of discrete Green functions, hence the determination of the time 
reversal of discrete random walks; then one takes the continuous limit. 

Proposition 2. The time reversal of the RBM with reflection angle —a starting 
from and stopped when hitting {z : 5Rz > x} is the RBM with reflection angle 
+a starting from the uniform distribution on the segment [x— i tana; x+i tana], 
conditioned not to hit {z : > x} again, and killed at 0. This last process is 
intertwined with a 3- dimensional Bessel process. 

Loosely speaking, the time reversal of a reflected Brownian motion in C a 
with reflection angle —a starting from is the reflected Brownian motion re- 
flected at angle +a, killed at 0, starting from "the uniform distribution at 
infinity" . 

Proof. Recall that we considered a random walk X on T intertwined with a 
random walk Y on N. Let G be the Green function for X killed when it hits 
{z : 3tz > TV cos a} for some large N. Let G be the Green function for Y killed 
at level N. Then it is obvious from the intertwining of X and Y that: 

G(0, x cos a + iy sina) = G(0, x) 

where < x < N, —x < y < x, x, y € N. As previously mentioned, the function 
i \— > — W is harmonic for the Markov kernel of Y (except in 0). From this, it is 
easy to compute: 

G{0,y) = 2 — 

cos z a 

Then Nagasawa's formula (sec e.g. RogWil93 , HI. 42) yields the Markov kernels 
of the time reversals X and Y of X and Y . For instance, we record the transition 
probability for X, if x G T is an inner point: 



p(x, x + i sin a) = p(x, x 
p(x, x — cos a) = 
p(x, x + cos a) = 



1 N-x 
sin 2 a N — x— 1 
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Notice that, as N goes to infinity, one gets the original transition probabilities. 
On the boundary, the reflection angle is reversed, which is not surprising. At 
this point, one takes the continuous limit as in section 1 (obviously the time 
reversal operation is compatible with the continuous limit) . The various claims 
follow easily. □ 

The conformal equivalence zhz^ maps C a onto Cp a . After an appropriate 
time change, this yields a general result for the time reversal of a reflected 
Brownian motion in an infinite wedge. 

In the set-up of the proposition, one may notice that up to its first hitting 
of the boundary, the time-reversed RBM has the law of a Brownian excursion 
in the half-plane {z : 5Rz < x} (for background on Brownian excursions, see 
|LSW02llV03| 1. One may use this to compute some probabilities of (indirect) 
interest in critical percolation. For simplicity, we will consider a RBM Z in the 
cone A = (1/2 + iy/3/2) — iC^fe starting from the apex and stopped on hitting 
the real line at time t (we are looking at a RBM in an equilateral triangle). We 
have seen that Z T is uniformly distributed on [0, 1]. Let g be the last time spent 
by Z on the boundary before r. We are interested in the joint law (Z g , Z T ), or 
rather in which side was last visited by Z conditionally on the exit point Z T ; so 
we will consider the event: 

K = {VIZg > 1/2} 

Since the time reversal of Z is a Brownian excursion until it hits the boundary, 
we have to compute the "harmonic measure" for the Brownian excursion. 

Let 7 be a Brownian excursion starting in H = {z : 3z > 0}, and let T be 
the first time it hits dA. Let e > 0, M > 0, and B a complex Brownian motion 
exiting the strip {z : < 3z < M} at time Tm- As M goes to infinity, the 
Brownian motion B conditioned on exiting the strip by the top converges to the 
Brownian excursion. Hence, making use of the Markov property of B: 

Px+fe(^T e dy) = lim P x+is (B T G dy\$B TM = M) 

M— >oo 



= F x+l£ (B T G dy) lim 



» v (%B Tm = M) 



m^oo F x+iE ($B TM = M) 
= ^-¥ x+lE (B T G dy) 

It is well known that harmonic measure is conformally invariant (see e.g. 
|RevYor!JT] ). and that the harmonic measure on the real line seen from a + ib is 
given by a Cauchy law (see e.g. |RogWil93| ) with density: 

bdx 

ir(b 2 + (x - a) 2 ) 

Moreover, according to the Schwarz-Christoffel formula (see |Ahl79| ) . the holo- 
morphic map: 

F(z) = I" u- 2 / 3 (l - u)- 2 ^ 3 du/B(l/3, 1/3) 



(i 
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is the conformal equivalence between the upper half-plane H and the equilateral 
triangle Afli that maps (0, 1, oo) to (0, 1, 1/2 + iy/3/2). If a £ (0, 1) is such 
that F(a) = x, we see that: 

Px(Y T e dy) = % vf-u \ ^2 

irF'(a){F 1 (y) - a) 2 

Then, one may compute: 

^ // x i z" 00 * /■'(«(«- !)r 2/3 V3 , 

P(K|z T e = F' (a) - ^ _ - ^ 4pL T i 

= F , (o) _x vs r *. 



2ttB(1/3, 1/3) 7i («(«- l)) 2 / 3 (u-a) 

7T\/3 



(a(l-a)) 2 / 3 2 Fi (1,4/3; 5/3; o) 



3r(2/3) 3 



B(2/3,2/3)7 (t(l-t))V3 

where 2-Pi(l, 4/3; 5/3; .) designates a generalized hypergeometric function (see 
eg |Bat53| or |AS65j . especially formulas 15.3.1 and 15.2.5). 



5 Relation with Watts' formula 

In this section, we recall Watts' formula ( |Wa96| ) and explain how it may be 
translated in the SLE language. This relation was suggested to us by Wendelin 
Werner. We shall see that it is closely related with the formula we derived 
above; in fact, this is one of the initial motivations for the present paper. 

Recall that Cardy's formula [022] gives the asymptotic behaviour of the 
probability of an open crossing between two sides of a topological rectangle in 
the limit when the mesh of the lattice goes to zero (for critical percolation) . This 
is equivalent to the fact that the hitting distribution of the exploration process 
in an equilateral triangle is uniform. This was shown by Smirnov |Smif)l| to 
hold in the case of site percolation on the triangular lattice, and it follows 
that the scaling limit of the whole exploration process is SLEg. Note that it is 
not difficult (see |LSW01p to prove directly in the continuous setting that this 
hitting distribution for SLEg is uniform. 

Trying to generalize Cardy's results using the same approach based on con- 
formal field theory, Watts |Wa96| considered the event that there exists simul- 
taneously an open left-to-right crossing and an open top-to-bottom crossing of 
a topological rectangle. He proposed a formula to describe the asymptotic be- 
haviour of the probability of this event when the mesh of the lattice vanishes, 
that seems to fit well with numerical simulations. Just as in the case of Cardy's 
formula, the double-crossing event can in fact be rephrased in terms of the 
exploration process. 
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For simplicity, we will discuss critical percolation on the triangular lattice (so 
each vertex is colored in black with probability p = 1/2). Let R be a topological 
rectangles, its boundary consisting of four disjoint arcs: 8l, &r, d R , d R (left, 
top, right, bottom). Let C(di,dj,c) be the event that there exists a crossing 
between di and dj with color c (here i, j G {L, T, R, B}, and the color c is black 
or white: c € {b, w}). We will also need the event T(9j,c) that there exits a 
connected component with color c linking the four boundary arcs except maybe 
di. As the triangular lattice is self-matching (see |Kes82) L it is classical that 
the two events C(8l, d R , w) and C(dr, d R ,b) are complementary: 

C{d L ,d Rl w) = C(<9 T ,c> B ,&) 

Now we are interested in C(3l, Or, 6)0(7(9^, 9b, b). With a little plane topology, 
one sees that: 



C(d L , Or, b) n C(dr, d Bl b) = T{d Tl b)\C{d Ll d R , w) 
T(dr,b) = C(d L ,d R ,b)\T(dr,w) 

Figure 4 illustrates these relations. 




Figure 4: Paths in a topological rectangle 
Translating this into probabilities: 

F(c(d L ,d R , b) n c(d T ,d B ,b)) = P(T(d T , b)) - P(T(cV, b) n c(d L ,d R , w)) 

and 

P(T(d T , b)) = F(C(d L ,d R , b)) - P(C(3 L , d R , b) n T{d T , w)) 

Keeping in mind that each vertex is colored with probability 1/2, so that chang- 
ing all the colors is measure-preserving, one gets: 

P(C(5 L , d R , b) n C(d R , d L ,b)) = P(C(d L ,d R , b)) - 2¥(C(d L , d R , w) T(d T , b)) 

Consider now the situation in an equilateral triangle ABC. Let X be the 
rightmost point on BC that is separated from AC by a black path. Cardy's 
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formula tells that X is uniformly distributed on BC (see figure) . The topological 
rectangle ADXE delimited by the upper- half of AB, the rightmost black path 
(solid), the leftmost white path (dashed), and the upper- half of AC is either 
crossed by a white path from AD to XE or by a black path from DX to EA. 
The exploration path goes from A to X leaving white vertices on its left and 
black vertices on its right. In the case where there is a white crossing between 
AD and XE, it is clear that the last edge of ABC visited by the exploration 
process before it reaches X is AB. Conversely, if there is a black crossing 
between DX and EA, the exploration process last visits AC. 



A 




Figure 5: The -stylized- exploration process (dotted) 

Smirnov has proved that the exploration process converges to a chordal SLEg 
(see ISmiOlp . So let 7 denote the trace of a chordal SLEg from A to B in ABC , 
t be the first time it hits BC, and g the last time it hits (AB) U (AC) before t. 
We may now state the SLE interpretation of Watts' formula in the topological 
rectangle ABXC (the left-hand should be understood as the scaling limit of the 
corresponding discrete probabilities): 

P(C(AB, XC, b)nC(BX, CA, b)) = P( 7r e BX)-2 f P(j g e AC\ 1t = Y)df( lT 

Jo 

Recall that ¥(j T e BX) = BX/BC (Cardy's formula). Watts' conjectural 
formula is equivalent to: 



^ &AC ^- X ^mkw)l 



BX/BC df 



(t(l-t))V3 
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Let us explicit the relation between chordal SLEg in ABC and the RBM we 
studied in the previous section. As mentioned in the introduction, the hull of 
chordal SLEg in ABC from A to B stopped at time r and the hull of the RBM 
in ABC stopped at its first hitting of BC are identical in law. Though, the 
processes producing these hulls are starkly different (see |LSW02p . Indeed, the 
SLE trace is non self-traversing, while the RBM is likely to go back through its 
past hull. Let Z be such a RBM, t the first time it hits BC, g the last time 
before r it visits (AB) n (AC). We have seen that: 

1 rBX/BC dt 
¥(Z a e AC\Z T =X) = — — / — — — - 

19 ' ' 5(2/3,2/3) J - t))V3 

and this is exactly the formula proposed by Watts for the scaling limit of critical 
percolation clusters, i.e. SLEg. Note also that the lowest points of the hulls 
on AB, say, have the same law since they are hull measurables. Since 7 is non 
self-traversing, j g is either the lowest point on AB or the lowest point on AC; 
but a.s. Z g is not the lowest point on AB or AC. 

Given these differences between the underlying curve laws, replacing the 
SLEg by the reflected Brownian motion in the "last visited edge" probability 
problem, the answer should a priori be different. Hence, there are two possibil- 
ities: 

• For some reason, the SLEg and RBM "last visited edge" probabilities 
coincide, and Watts' prediction holds. 

• Watts' prediction does not hold. 

6 Vases 

In this section, we generalize the previous properties of RBM, replacing wedges 
by more general domains with a symmetry axis. More precisely, let h : R + — > R + 
be a positive, diffcrcntiablc function. Suppose h(0) = 0, and h[x) > for x > 0. 
Then we will consider the "vase" Ch- 

C h = {zeC : diz > 0, |9te| < h(U(z))} 

The shape function h is fixed in this section, so we may omit the subscript. 
As in the previous section, we define a tailor-made discrete model to prove an 
intertwining relation for its continuous limit. The situation here is slightly more 
complicated, so we will use a continuous time, discrete state space Markov chain. 
For the sake of simplicity, suppose the function h is strictly increasing. Let N 
be some large, fixed integer, and let Xk = h^ 1 (k/N). We will consider a graph 
r (depending on h, N) with vertices [xu +iy 7 / N)k^n,ye , z,\y\<k- For k € N, define 
cik = arctan(l/(a;fe+i — Xk))- Let us now define the jump rates (non diagonal 
elements of the Q-matrix). If z = Xk + iy/N is an inner point of T (i.e \y\ < k), 
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the jump rates are : 

Q(z, z + i) = Q(z, z — i) = \ 

Q(z,Xk+i +iy/N) = [cotanafc(cotanafc + cotana^^!)] -1 
Q(z, Xk-i +iy/N) = [cotanafe_i(cotano!fc + cotanafc_i)] _1 

These jump rates are chosen in order to get a zero mean, isotropic walk. For a 
boundary point z = Xk + ik/N: 

Q(z,z-i/N) =1/2 
Q(z, x k+1 + ik/N) = Q{z, x k+1 + i(k + 1)/N) 

= [cot an (cot an a & + cotanafc-i)] -1 

For the apex, set Q(0,xi) = Q(0,xi ± 1/iV) > 0. Here the jump rates are 
chosen so that intertwining relation holds. More precisely, let A be the Markov 
transition kernel from (x k )kefi to T defined as: 



a(^.) = 2^t£> 



k + iy/N 

y=-k 

Then the intertwining relation is valid for Q-matriccs. More precisely, let Q be 
the Q-matrix of the Markov process just defined, and Q be the Q-matrix of the 
Markov process on {.T^jfegN C K.+ defined by the jump rates: 

{Q(xk, Xk+i) = Hiy[cotanafe(cotanai, + cotano^i)] -1 
Q(xk, Xk~i) = 2k+i [ c °tauQfc-i(cotanafc + cotanafc_i)] _1 

for k > 0, and Q(0, xi) = 3Q(0, Xi). It is immediate to check that: 

AQ = QA 

Then it is sufficient to exponentiate this intertwining relation (in appropriate 
Banach spaces, say End(7i(r)) and End(Zi({a:fc}fceJv)) ) to get the intertwining 
relation for the associated semigroups (-Ft), (-Ft): 

AP t = P t A 

for all positive time t. We now examine the scaling limit of the second process. 
Let / be some arbitrary function in C? (R+). Then, as N goes to infinity, if 
k = kpf is such that x kN — ► x for some x € R + : 

1 / 2k + 3 

Qf(x k ) = , o; ; ; r — (/(Zfc+l) - f{x k )) + 

(2k + l)(cotanafc + cotanafc_i) \cotanafe 

9k — 1 

-(f(x k ^)-f(x k )) 



cotanafc- 



15 



After a time rescaling, these discrete generators "converge" to the diffusion 
generator: 



ti d Id 2 
h dx 2 dx 2 

One checks immediately that a scale function for this diffusion on R + is: 



For a wedge, h{x) = xtana, then x > a;" 1 is a scale function for the associated 
diffusion, i.e. the 3-dimensional Bessel process. 

Wc now discuss the proof for the convergence of the processes defined above. 
It is similar to the previous proof, with added technicalities. Firstly, one has 
to adapt the submartingalc-problcm characterization for vases. One possibil- 
ity is to map conformally Ch to a wedge, say C v /4 (for the Riemann Mapping 
Theorem, see |Ahl79| ): after a time change (see |RevYor9"T] for a discussion of 
the image of a complex Brownian motion under a conformal map), one may 
use a variant of the Varadhan- Williams result. Here the reflection angle may 
vary along the boundary, but stay negative. Then one argues as above: sub- 
sequential scaling limits exist (tightness); they satisfy a submartingale-problem 
-a martingale problem for the one-dimensional diffusion-, as seen by dominated 
convergence after applying Dynkin's formula; hence the limit is uniquely de- 
termined, and the processes converge weakly. Reflections are taken care of as 
before, via an occupation time estimate. The intertwining is obviously preserved 
in the continuous limit. Hence one may state: 

Proposition 3. The reflected Brownian motion in the vase Ch with constant 
reflection index —1 is intertwined with the diffusion on K + corresponding to the 
generator \-§^ + \ -§^i by the Markov transition kernel from K+ to Ch: 

A(x, .) = Unif ([x — ih(x);x + ih(x)\) 

Let us stress that the reflection angle varies along the boundary; the reflec- 
tion at x + ih{x) points in the — i exp(2z arctan/i'(x)) direction. Starting from 
an arbitrary diffusion on E + , one may consider an increasing scale function </>. 
Setting h(x) = l/^(f>'(x), one constructs an intertwined reflected Brownian mo- 
tion. Of course, one can also consider shape functions defined and positive on ]R 
("funnels" rather than vases), or on a segment [a, b], with h(x) = if and only 
if x E {a,b} (closed vases). We will now consider an example of the previous 
construction. 

Corollary 1. Consider the shape function h(x) = x 13 , with f3 > 0. Then the 
real part of the RBM in the vase Ch with constant reflexion index —1 is a Bessel 
process of dimension 2(3 + 1. 

When the shape function h is not monotone, one has to be a bit more 
cautious about the definition of the grids, but this requires no substantial 
changes in the proof (the terms 2k + 1 are to be replaced by l(xk), with 
l(x k+ i) - l(x k ) = 2sgnh'(x k ))- 
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